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CORONAE AS CONSEQUENCE OF LARGE SCALE MAGNETIC FIELDS IN TURBULENT ACCRETION DISKS 
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ABSTRACT 

Non-thermal X-ray emission in compact accretion engines can be interpreted to result from magnetic dissi- 
pation in an optically thin magnetized corona above an optically thick accretion disk. If coronal magnetic field 
originates in the disk and the disk is turbulent, then only magnetic structures large enough for their turbulent 
shredding time to exceed their buoyant rise time survive the journey to the corona. We use this concept and a 
physical model to constrain the minimum fraction of magnetic energy above the critical scale for buoyancy as 
a function of the observed coronal to bolometric emission. Our results suggest that a significant fraction of the 
magnetic energy in accretion disks resides in large scale fields, which in turn provides circumstantial evidence 
for significant non-local transport phenomena and the need for large scale magnetic field generation. For the 
example of Seyfert AGN, for which ~ 30 per cent of the bolometric flux is in the X-ray band, we find that more 
than 20 per cent of the magnetic energy must be of large enough scale to rise and dissipate in the corona. 

Subject headings: accretion, accretion disks — black hole physics — instabilities — MHD — turbulence 



1. INTRODUCTION 

Accretion disks are widely appreciated to be a powerful 
source of emission from gas or plasma orbiting central stars 
or compact objects (see, e.g.. iFrank et al.ll2002f) . In order to 
explain the rapid variability and short lifetimes of accreting 
systems without unphysical mass densities, some enhanced 
angular momentum transport beyond that which can be sup- 
plied by the microphysical transpor t coefficients is typically 
required (Shakura & Sun yaevl[l973[) . Many accreting sources 
show jets, outflows, and active coronae highlighting that disk 
dynamics and energy release involves some combination of 
local and large scale transport. Understanding the relative bal- 
ance between the two is of fundamental importance. 

In some disk models, the primary angular momentum trans- 
port and dissipation takes place above the disk (Lynden-Bell 
1969; Field and Rogers 1993) and the turbulence within the 
disk itself plays a secondary role. However, the magnetoro- 
tational instability (MRI) has emerged as a likely source of 
turbulence within accretion disks, and a leading candidate to 
contribute local turbule nt angular momentum transport in suf- 
ficiently ionized d isks ( Velikhovl l 19591 : fchandrasekhari l 1 9601 
Balb us & Hawlevlll99lL 119981). Th ree-di mensional numeri- 
cal simulations |Hawley et alj|1995l, 1 1996t [Brandenburg et all 
ll995tlStone et al.ll 19961) have revealed that the nonlinear evo- 
lution of systems unstable to the MRI leads to sustained MHD 
turbulence and outward angular momentum transport. 

Understanding the saturatio n of the MRI is curren tly a 
topic of active research (see iPessah & Goodman! 120091 and 
references therein). Local shearing box simulations have 
not converged to practical angular momentum transport co- 
efficients, which are found to depend strongly on the sim- 
ulat ion box size and the initial seed magnetic field (see , 
e.g.lHawlevetalJll995l: ISano et al.ll2004t IPessah et all 120071: 
iBodo et all 120081), as well as the mag netic Prandtl number 
(Flem ing et al.1 120001: ISano & I nutsuka 20011: iFromang etail 
l2007tlLesur & Lon garetti 12007b . An important frontier in this 
regard is to understand the development and role of large scale 
magneti c field structures in the saturation of MRI-drive n tur- 
bulence dMiller & Stond l2000; Bra ndenburg et al.ll 19951) and, 



more generally, in the process of angular momentum trans- 
port. A realistic disk like ly involves coupled inter nal and 
coronal dynamics (see, e.g jKuncic & Bickne ll 2004). 

Observationally, the relevance of large scale magnetic fields 
in accretion disks is strongly motivated by the interpretation 
of X-ray flux in Seyferts which has been best interpreted as 
coronal emission. The flux from 1-500 keV ranges from 10 
to 50 per cent of the total flux (Mushotzky et al. 1993). 
Galactic black hole X-ray sources show both thermal and non- 
thermal (power law) spectral components, with the ratio of 
non-thermal to total luminosity ranging between 20 and 40 
per cent (Nowak 1995). The leading paradigm for X-ray emis- 
sion in these accreting systems involves an optically thin, hot 
corona powered by magnetic field dissipation (e.g. Haardt & 
Maraschi 1993: Field & Rogers 1993). If the corona results 
from magnetic structures dissipating above the disk midplane 
that were originally produced within the turbulent disk (e.g., 
via magnetic instabilities, such as the MRI), then these struc- 
tures must be of large enough scale to survive the buoyant rise 
without being prematurely shredded by disk turbulence. 

If all coronal and jet emission results from fields initially 
produced within a turbulent disk, then the fraction of coro- 
nal to bolometric luminosity is directly related to the fraction 
of magnetic energy associated with buoyant fields of large 
enough scale to survive the vertical trip without being tur- 
bulently shredded. In this Letter we employ this concept to 
develop a model relating the observed ratio of coronal to bolo- 
metric emission to the fraction of magnetic energy produced 
in the disk that is of large enough scale to buoyantly rise to 
the corona. By comparing the implications of our model with 
observations we infer that a significant fraction of magnetic 
energy in accretion disks must reside in large scale fields. 

2. WHY CORONAE REQUIRE LARGE SCALE FIELDS 

In order for magnetic fields to power coronae and jets, the 
buoyancy time, associated with a magnetic structure rising 
through the disk must be smaller than the time associated with 
its turbulent diffusion, id- These timescales are estimated as 
th = H/Ub, where H is the disk half thickness and Ub is 
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the characteristic buoyancy speed, and t& = I jvt , where 
vt ~ vlt is the turbulent magnetic diffusion coefficient, v is 
the dominant turbulent speed, and l t is the characteristic scale 
of a typical (anisotropic) turbulent cell, i.e., l t ~ (Ixlylz) 1 ^ 3 - 
The escape condition, t& < sets a lower bound on the 
scale I of magnetic structures that can survive shredding and 
reach the coronae (Blackman & Tan 2003), namely 

v 

Ub 



I 2 >l 2 = l t H-^r- 



(1) 



The critical scale applies to the smallest dimension of a 
given magnetic structure, i.e., l c < mm{l x ,l y ,l z }, since 
this determines the shortest diffusion time. For a thin mag- 
netic loop / would be the loop thickness not the distance 

between foot-points. In the a-viscosity disk framework 

dShakura & Sunvaevl [19731) . v t = vl t = ac s H ~ v 2 /n, 
where c s is the sound speed and fi is the local angular fre- 
quency. Using c s ~ flH, this implies v = a 1 ^ 2 ^ and thus 
It ~ a 1 / 2 !!. Plugging these into Equation (Q} implies 
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In a fully developed MHD/MRI turbulent flow the velocity 
of energy containing turbulent motions is approximately (if 
not slightly less than) the RMS Alfven speed, i.e., v ~ va- 
Therefore if U b < v A ~ v ~ a 1/2 c s , then I > l c > o^^H = 
It /a 1 / 4 , which is larger than It for a < 1. If instead we 
use Ub(l) < Cs, the analogous procedure gives I > l c > 
a 1 / 2 !! = l t . Either way, the buoyant magnetic structures that 
survive turbulent shredding in the disk must have I > l c > l t . 

That l c equals or exceeds the characteristic turbulent scale 
implies that an accretion engine with a significant fraction of 
power emanating from coronae requires a significant fraction 
of magnetic energy to be organized in magnetic structures of 
large scale with I > l c . In the sections that follow, we quantify 
this fraction by determining l c and Ub(l) and connect them 
with the observed ratio of coronal to bolometric emission flux. 

3. CORONAL EMISSION FRACTION 

We develop a model in which the observed fraction of coro- 
nal to bolometric disk luminosity is determined by the rate of 
large scale magnetic energy rising to the corona. We assume 
that buoyant structures fill a volume fraction f v in the disk 
such that the average disk mass density is given by 



P = Po - fv(Po - Pi) ■ 



(3) 



where p a and pi are the mass densities external to and inter- 
nal to buoyant structures, respectively. The average magnetic 
energy density is then 
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(B 2 Q + f v B 2 s ) 



(4) 



where B 2 S = B 2 — B 2 corresponds to the difference of in- 
ternal and external magnetic fields squared. We suppose that 
although the entire magnetic energy of the structures Bf con- 
tributes to their initial buoyancy, only B 2 S survives to the 
corona; the smaller scale fields are "bled" away during the 
buoyant rise. We define the fraction of the magnetic energy 
density in scales larger than the critical scale for surviving the 
buoyant rise to be 
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The factor f v arises in the numerator because B 2 S is contained 
only in the volume of the buoyant structures. The quantity / s 
represents the fraction of magnetic energy with scales larger 
than l c , which can be written more generally in terms of inte- 
grated magnetic spectra as 



f^ in E M (k)dk 
J^ x E M (k)dk 



(6) 



where the limiting wavenumbers are k c — 2ir/l c , fc m j n = 
2tt/H, and fc max = 2-7r/7diss and Zdiss is the dissipation scale. 

Motivated by the physical picture described above, we 
break up the accretion energy per unit area dissipated at a 
given radius into the sum of the dissipation associated with 
small scale field within the disk 



D d = QY,n 2 v t {l ~ fs) = 2Qa/3c s e mag (l - /. 
and the dissipation of large scale field in the corona 



(7) 



(8) 



Here, S is the surface density, Q = (d\nil/dR) 2 /2 and 
j3 = pc 2 /2e mag , where e mag = B 2 /8ir is the total magnetic 
energy density. The first term on the right of Equation (|7]i re- 
sembles that whic h would follow from standard disk theory 
dFrank et al.ll2002l) but with the extra factor of 1 - / s . 
Defining the ratio of coronal to total dissipation as 

the ratio of coronal to disk dissipation becomes 

IX _J_ fs Ub 

D d 1-q 2Qaf3(l - /„) c s ' 

and depends cleanly on the ratio U b /c s . Since both D c and 
D d are expected to be dominated by their contributions near 
the inner radius, we do not address the radial dependence of 
these quantities in detail here and interpret q as an estimate of 
total coronal to bolometric emission. 

4. CONSTRAINING THE BUOYANCY SPEED 

To constrain the ratio Uh/c s , we consider forces on a mag- 
netic structure in pressure balance with its exterior, i.e., 



(10) 



B 2 B 2 

+ rnk^T, = —2- + n k B T . 
sir sir 



(11) 



where n and T are the corresponding number densities and 
temperatures, and ks is the Boltzmann constant. Assuming 
that Ti = T = T, we obtain 



Po - Pi 



87TC? 



(12) 



We take the force density acting on the magnetic structure 
to be F b — where F b is the upward gravitational buoyancy 
force density and F dl is the drag force density. The former can 
be estimated according to 



*b = gH- 



B 2 



8irH 



(13) 



where the last equality follows from hydrostatic equilibrium 
for an isothermal gas (T — 1), that is kBT/m p = c 2 = 
(GM/R)(H/R), with g = GM/R 2 . We approximate F d 
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FIG. 1. — Ratio of the critical scal e fo r buoyancy to the disk scale 
height, l c /H, that satisfies Eqs. J5) and 115\ simultaneously, as a function 
of the density ratio pi/po- The three sets of solid curves correspond to 
a = {0.5,0.1,0.01} from top to bottom. Larger values of a correspond 
to more efficient shredding and thus require larger scales l c to survive the 
buoyant rise. For a given a, the solutions correspond to a range of drag coef- 
ficients 10 — 3 < C(j r < 2, logarithmically spaced. The bottommost curves 
in each set correspond to the lowest drag; larger C^r requires a larger Z c to 
survive buoyant rise. The dashed lines show the turbulent scale It = a^^H. 
That each set of curves lies above the line associated with It (ce) shows the 
importance of large scale fields. A corona cannot form for region l c > H . 

to be the high Reynolds number hydrodynamic drag associ- 
ated to a cylinder of length L and diameter I, namely (see, 
e.g. Moreno-Insertis 1986; Landau and Lifshitz 1987) 
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with drag coefficient Cd r of order unity. The work done per 
unit volume by the net force over a distance H equals the 
kinetic energy density of the rising structure, H{Faj. — Fb) ~ 
PiU^/2. Combining this with Equations (fT~2l> — dl4T> we obtain 



= 1 



Pi 

Po 



1/2 (2C dx H 
\ ir I 



Pi 
2p 
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(15) 



Therefore, for all buoyant structures Ub(l) > U(l c ) since 
I > l c . Given that I < H for all structures that fit in the disk, 
Ub(l) < c s for all I would imply pi/p > 2/3 — 4Cd r /37r. 
We thus restrict ourselves to this density ratio regime. Note 
that our estimate for Ub/c s is itself an upper limit since 
we consider only a hydrodynamic drag force restricting the 
buoyant rise, ignoring for exampl e, magnetic tension, c.f. 
ISchramkowski & Torkelssonl d2002l) . Whether the dynamics 
allows densities below the above upper limit remains an open 
question. However, even if the structures were initially able 
to move faster than c s , we would expect shocks and the asso- 
ciated dissipation to slow the motion via additional drag. 

Vishniac (1995) estimated a buoyant rise time {/ b ~ ac s . 
This value is low compared to the upper limit c s we consid- 
ered. Had we used Ub ~ ac s , then the value of / s in Equa- 
tion (TTTb would dramatically increase; highlighting the im- 
portance of large scale field for coronal dynamics even more. 
Some characteristic values of / s resulting from this smaller 
buoyancy speed are shown as dashed lines in Figure [2] 

5. QUANTIFYING THE IMPORTANCE OF LARGE SCALE FIELDS 



FIG. 2. — Lower limits on the fraction of magnetic energy residing at scales 
large enough for their buoyant rise time to be small compared to the corre- 
sponding turbulent shredding time, / s , as a function of the ratio of coron al to 
bolometric dissipation, q. The solid curves correspond to U\, = c s in Eq. )16t 
which gives the most stringent lower limit for a/3 = {1.0, 0.5, 0.1}, from 
top to bottom. The dashed curves correspond to = a c s , as suggested in 
Vishniac (1995). In this case, the lower limit given by Eq. {16} is independent 
of a. The dashed lines correspond to j3 = {100, 10, 5}, from top to bottom 
((3 = 1 would coincide with the curve for a/3 = 1 for = c s ). This last set 
of curves show that the same ratio of coronal to bolometric emission requires 
a higher fraction of large scale fields as the buoyancy speed is reduced. 

Setting I — l c in Equation ( fTBT ), and using Equation (f2j), 
gives a 5 th order equation for l c /H as a function of Cdr, 
Pi/ Po, and a. The physical solution to this equation is shown 
in Figure[T]for a range in drag coefficients, 10~ 3 < Cd r < 2, 
and for three values of a = {0.5,0.1,0.01} as a function 
of the density ratio p\j p - The characteristic turbulent scale 
l t = a 1 / 2 !! is also shown as horizontal lines for each value 
of a. Since l c is the minimum scale for buoyant rise, the fact 
that l c > It highlights the importance of large scale fields. For 
the regime on the plot where l c > H, a buoyant corona cannot 
arise from fields produced internally to the disk. 

We can set two constraints on the minimum fraction / s of 
magnetic energy that the disk must produce in fields with 
scales I > l c for a given coronal to bolometric emission 
fraction q. The more stringent bound relies on the fact that 
Ub{l) < Ub{H) for / < H; a less severe limit is obtained 
requiring Ub(H) < c s , which is satisfied as long as the den- 
sity ratio pi I p does not fall below the lower limit discussed 
above. Applying these conditions to Equation ( TTOb we obtain 
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f s U h (H)/c s 
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(16) 



q 2Qa/3(l - /„) " 2Qa/3(l - /.) 

Then, for an observationally inferred value of coronal to bolo- 
metric flux, we can obtain a lower limit on the fraction of 
magnetic energy residing in large scale fields, i.e., 

a(3 



fs> 



a/3 + 4(l/?-l)/9 ' 



(17) 



where we have used that Q ~ 9/8 for a Keplerian disks. 

The lower limit for / s depends on the dimensionless pa- 
rameters characterizing the angular momentum transport effi- 
ciency, a, and magnetic pressure support, (3, only though the 
product a(3. This result is very encouraging because, despite 
the fact that both quantities vary over several orders of mag- 
nitude across simulations carried out in domains with various 
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sizes and with diff erent field strengths and geometries (see 
Pessah et al. 2006a, and references therein), their product re- 
mains nearly constant 1 with a/3 ~ 0.5 (see Blackm ail et al.l 
2008, and references therein). 

The upper limits for the fraction of magnetic energy associ- 
ated with large scale field structures are shown in Figure|2]for 
three values of the product a(3 — {1.0,0.5,0.1}. The min- 
imum constraints on / s illustrated show that, if the observed 
non-thermal emission is interpreted as coronal emission due 
to magnetic dissipation from buoyant fields that were pro- 
duced within a turbulent disk, then a significant fraction of 
the energy budget of the magnetic field build in the disk must 
be produced in fields of scale I > l c . Since Figure Q] shows 
that in general I > l t , together these figures highlight the im- 
portance of large scale magnetic fields in powering coronae. 

Finally we note that MHD j e ts mo dels typically invoke 
global scale fields, c.f. iFerreiral ([2007). If these fields arise 
from the opening of coronal fields (as in the sun; Wang & 
Sheeley 2003: Blackman & Tan 2004) then the mechanical 
luminosities of jets would represent an additional contribu- 
tion to that which results from the buoyant rise of magnetic 
fields. Specifically, D c in Equation (0 would be repalced by 
D c + Dj where the latter is the jet power. This would further 
increase our lower limits on f s . 

6. RELATION TO PREVIOUS WORK 

In our physical picture, the coronal emission fraction q de- 
pends on the fraction of magnetic energy / s produced in scales 
larger than the critical scale l c . We incorporate the density 
contrast between buoyant structures and the ambient medium 
required for coronal feeding. Also, our lower limit on / s em- 
ploys c s as an upper limit for the buoyant rise time. We use the 
a-viscosity prescription only for dissipation inside the disk; 
the coronal dissipation is modeled as a distinct contribution 
(see Eqs. Q and J8|). These features differ from those of 
Merloni & Fabian (2002) and the ones summarized in Wang 
et al. (2004). In those treatments, there is no explicit dis- 
tinction between the density inside and outside the buoyant 
structures or the role of large and small scale magnetic fields. 
Furthermore, the coronal emission fraction is considered to 
be a subset of the total dissipation, modeled entirely with the 
a-viscosity prescription. Also, the buoyant rise time is taken 
to be the Alfven speed; less than our upper limit value of c s . 

In the present work we implicitly consider systems with low 
enough accretion rates such that radiation pressure is unim- 
portant. A subtlety associated with radiation pressure is that 
the thermal photosphere can be significantly higher than the 

1 The constancy of a/3 is consistent with the relations below Equation (2); 
for a disk with MRI growth time ~ Q^ 1 we expect ut = cec B H ~ v 2 /Q ~ 
vj^/fl. Using c s ~ Q.H then gives a oc /3~ 1 with a proportionality constant 
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scale at which the magnetic pressure dominates the thermal 
pressure (e.g. Hirose et al. 2009). Thus the non-thermal coro- 
nal emission arises at larger scale heights compared to when 
radiation pressure is ignorable. The reduction of the coronal 
emission fraction for such large accretion rate systems is ob- 
served in Wang et al. (2004). For a fixed magnetic spectrum, 
this would also be expected in our paradigm, because H in 
Equation (O is the scale that we consider a buoyant struc- 
ture must rise to contribute to coronal emission. If the buoy- 
ant structure has to move higher, then l c would be larger and 
less of the magnetic energy would survive the buoyant rise. 
In Merloni & Fabian (2002) the reduction in coronal emis- 
sion for large radiation pressure occurs because their coronal 
emission fraction depends inversely on the total disk pressure. 

7. CONCLUSIONS 

Starting with the assumption that coronal luminosity from a 
turbulent accretion disk results from buoyant magnetic struc- 
tures that survive turbulent shredding for at least one vertical 
density scale height, we derived lower limits on: (i) the scale 
of such magnetic structures and (ii) the fraction of magnetic 
energy that needs to be produced above this scale within the 
disk to account for observed values of coronal to bolometric 
luminosity. In our minimalist model, we considered the buoy- 
ant structures to be in pressure equilibrium with the ambient 
medium but to have an additional magnetic energy contribu- 
tion from scales above the critical scale l c , and a lower density. 

We find that typical ratios of coronal to bolometric luminos- 
ity observed in AGN require the critical scale for buoyancy to 
robustly exceed the characteristic scale set by turbulent mo- 
tions and that double digit percentages of magnetic energy 
should reside in fields above this scale. This is consistent 
with recent work highlighting the importance of in situ large 
scale dynamos in feeding coronae (Blackman 2007; Vishniac 
2009). Our results complement growing motivation to con- 
sider larger domains in stratified MRI simulations and moti- 
vate analysis of the magnetic energy spectra produced therein. 
The results also resonate with models of accretion disks in 
which buoyancy and coronal dissipation play a primary role 
for transport (Lynden-Bell 1969; Field & Rogers 1993). 
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